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ABSTRACT 

Using a sample of 100 HI- selected damped Ly-a (DLA) systems, observed with the High 
Resolution Echelle Spectrometer on the Keck I telescope, we present evidence that the scatter in the 
well-studied correlation between the redshift and metallicity of a DLA is largely due to the existence 
of a mass-metallicity relationship at each redshift. To describe the fundamental relations that exist 
between redshift, metallicity and mass, we use a fundamental plane description, which is described by 
the following equation: [M/H] = (-1.9±0.5) + (0.74±0.21)-logAv 90 - (0.32±0.06)-z. Here, we assert 
that the velocity width, Avg , which is defined as the velocity interval containing 90% of the inte- 
grated optical depth, traces the mass of the underlying dark matter halo. This description provides 
two significant improvements over the individual descriptions of the mass-metallicity correlation and 
metallicity-redshift correlation. Firstly, the fundamental equation reduces the scatter around both 
relationships by about 20%, providing a more stringent constraint on numerical simulations modeling 
DLAs. Secondly, it confirms that the dark matter halos that host DLAs satisfy a mass-metallicity 
relationship at each redshift between redshifts 2 through 5. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: fundamental parameters — 
galaxies: ISM — galaxies: kinematics and dynamics — quasars: absorption lines 



1. INTRODUCTION 

The study of high-redshift quasars has resulted in 
the detection of over 6000 damped Lyman a absorp- 
tion systems (DLAs), which are absorption systems 
with a n H I column densi t y greater than 2 x 10 20 
cm~ 2 (IProchaska et all l20l)5l IProchaska k Wolfd l200i 
iNoterdaeme et al.l 120091 I2012D . DLAs are thought to 
be th e progenitors of modern day galaxies ([Wolfe et al.l 
119951 ) , and therefore provide an important observational 
constraint on the formation and evolution of galaxies 
seen today. Most DLAs are found using low resolu- 
tion ins truments such as the Sloan Digital Sky Survey 
(SDSS; lAbazaiian et al.l 120091) . Although these instru- 
ments allow for accurately determining the column den- 
sity of the neutral hydrogen, follow-up measurements of 
DLAs using high resolution instruments are needed to 
accurately measure the resonance line transitions of met- 
als that are present in the DLAs. These measurements 
yield the column densities of the DLA's gas -phase con- 
stituents (for a review see: IWolfe et all 120051) . and allow 
for the determination o f the metallicity of the gas (e.g. 
Prochask a et al.l 1200 3^ . the cooling rate of the neutra l 
gas due to [C II] 158 /xm e mission (4: IWolfe et al.l l2003l). 
and the dust-to-gas ratio ([Pettini et all 119941 ) . Besides 
these quantities, the velocity profiles of the metal lines 
also gives us an insight into the kinematics of the gas (e.g. 
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IProchaska k WolfdH997t iLedoux et al.lll998D . Together 
with the redshift of the DLA, these parameters individ- 
ually provide us insight into the composition and nature 
of DLAs through the distribution that each of these pa- 
rameters assumes. In turn, these data offer insight into 
the properties of gas-rich and presumably starforming 
galaxies at high redshift. Furthermore, we may obtain 
additional information by considering how these param- 
eters correlate to each other. 

There are two such correlations for DLAs that have 
been extensively examined. The first of these corre- 
lations is the evolution of metallicity with redshift; as 
redshi f t decreases, met allicity increases (IProchaska et al.l 
[20031: IKulkarni et all 120013 IRafelski et all 120121) . The 
second is the correlation bet ween the kinematics and 
the metallicity of t he DLA (IWolfe k Prochaskal 119981 : 
ILedoux et all 120061: IProchaska et al.l I2008D . This sec- 
ond relationship is thought to be caused by the un- 
derlying mass-metallic i ty relationship seen i n galaxies 
(iTremonti et al.l 12004 iSavaglio et all 120051: lErb et all 
2006: IMaiolino et all l2008) , where the kinematics of the 
DLA is assumed to be related to the virial speed, hence 
the mass of t he dark ma tter halo that hosts the DLA. 
iFvnbo et all (|2010t 1201 ID used this correlation to image 
DLAs in emission. By selecting higher metallicity sys- 
tems, they improved their rate of detection, which they 
assert is due to the larger mass of the DLA. 

One of the main features of both correlations is the 
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large intrinsic scatter, which is significantly larger than 
the observational uncertainty in the measurements. One 
possible cause for the large intrinsic scatter around either 
correlation is that the two correlations are dependent on 
each other, since they both correlate metallicity to a sec- 
ond parameter. To quantify the potential dependency 
of the two correlations, we can see if a linear combina- 
tion of the two correlations is able to reduce the intrinsic 
scatter. Such a linear combination is called a 'funda- 
mental' plane, because it traces out a plane in the three 
dimensional parameter space it describes. 

The occurrence of fundamental planes in astronomy 
is relatively common. Indeed, in recent years there 
have been several publications on the existence of fun- 
damental planes in a variety of astrono mical areas, from 
black hole a ctivity (iMerloni et al.1120031 ) , to that of radio 
magnetars (|Rea et al.l 12012'). To illustrate the advan- 
tage of using fundamental planes, we will consider one 
of the first fundamental planes, the fundamental plane 
of elliptical galaxies. Elliptical galaxies show a great 
diversity in their properties such as surface brightness 
and velocity dispersion. However, most of these prop- 
erties a re correlated to the luminosity of the elliptical 
galaxy (jKormendvl 119821 ) . Indeed, both velocity disper- 
sion and surface brightness show a correlation with lumi- 
nosity. However, the scatter around both correlations is 
larger than the obse rvational uncertainty . In an effort to 
reduce this scatter, IDiorgovski fc Davis! (|1987t) noticed 
that when the luminosity, velocity dispersion and sur- 
face brightness are plotted on a three-dimensional plot, 
the points line up on a plane with a scatter around the 
plane that is within the observational uncertainty. This 
planar description therefore showed that the majority of 
the scatter in either of the correlations was caused by the 
other correlation. As such, the plane description provides 
a more stringent constraint on theoretical models than 
each of the individual correlations. In addition, the co- 
efficients that describe this plane provide us insight into 
the galaxy population as a whole; deviations from the 
expected values could indicate variations in the initial 
mass function, dark matter fraction and s tellar popula- 
tion effects (see e.g. lGraves fc F aber 20l3). 

In this paper, we explore the possibility of combining 
the redshift-metallicity and mass-metallicity correlation 
into a 'fundamental' plane for DLAs. We find that indeed 
such a plane exists, and applying this plane decreases the 
intrinsic scatter seen in both correlations significantly. 
We discuss new insights provided by this fundamental 
plane and we explore the physical origin of such a plane. 

We conclude the paper by looking for the existence 
of this fundamental plane in current numerical simula- 
tions. Since most numerical simulations do not accu- 
rately trace all parameters involved; direct evidence of 
the fundamental plane is difficult to find. Instead we 
focus on projections of the fundamental plane; i.e. we 
primarily compare two parameter correlations and distri- 
butions of a single parameter. One distribution in partic- 
ular that we examine is the velocity width distribution. 
This distribution has two notable features, a very large 
median and a tail toward large velocity widths. Pre- 
vious modeling showed that only two models were able 
to recreate these two features of the distribution. The 
first mod el posed that DLAs are thick rapidly rotat- 
ing disks ([Prochaska fc WoIfelll997l) . whereas the second 



model posed that DLAs are disk progenitors in a stan - 
dard cold dark matter cosmology (|Haehnelt et al.llT998|) . 
The latter, however, erroneously excluded syst ems with 
low velocity widths (Prochaska & Wolfe 120101 ) . More- 
over, later numerical simulations based on ACDM models 
that i ncluded a more realistic transport of ionizing radi- 
ation (Ra zoumov et al.l 120061: iPontzen et all 12008) were 
unable to recreate the distribution of observed velocity 
widths. This inability of simulations based on ACDM 
models to produce the observed velocity width distribu- 
tion is attributed to their difficulty in producing enough 
DLAs that reside in d ark matter halos with large masses 
(|Pontzen et al.ll2008D . We explore this hypothesis with 
the current data set. 

This paper is organized as follows. In Section [2] we de- 
scribe our observational strategy and the reduction pro- 
cess we used. We explain the method used to measure 
the parameters described for each DLA, consider any bi- 
ases in the sample, and explore their effects on the dis- 
tribution of the parameters. Since a fundamental plane 
is dependent on the correlations between the parame- 
ters involved, it is crucial to explore these correlations. 
This is done in Section [3] In Section |4] we introduce 
the fundamental plane. Finally, we will discuss these re- 
sults and compare them to numerical models in Section[5j 
The sample discussed in this paper is the largest sample 
of DLAs for which the metal lines have been observed 
with a single high resolution ( R > 40000) instrument 
(Keck /HIRES : lVogt et aT1ll994h . and to our knowledge, 
this is the first paper that considers multi-parameter cor- 
relations, such as fundamental planes, for DLAs. 

2. SAMPLE SELECTION, PARAMETER 
DEFINITION AND SAMPLE BIASES 

In this section we discuss the reduction process used 
for the sample. We explain our sample selection and any 
biases that could influence our results. 

2.1. Observations and Reduction 

All of the spectra in this sample were obtained 
with the High Re solution Echelle Spectrometer (HIRES: 
iVogt et al.l 119941 ) on the Keck I 10m telescope over the 
course of almost two decades of observing. Table Q] shows 
a detailed journal of observations. All observations prior 
to 2004 were carried out with the Tektronix 2048x2048 
CCD; the remaining observations were made with a 3- 
chip mosaic of MIT-LL 2048x4096 CCDs. The FWHM 
resolution of each object, Av res , varies depending on the 
specified slit size and the atmospheric conditions, but 
most of the data were taken with a 0.86" or 1.15" slit, 
which results in a maximum Av rcs corresponding to 6 
km s _1 , and 8 km s _1 respectively. Wavelength cover- 
age was dependent on the redshift of the DLA. Finally, 
we computed the average signal to noise ratio (S/N) per 
1.4 km s _1 pixel for each observation, by taking the cen- 
tral 200A of each observation and finding the median 
S/N for this range. If the range included too many ab- 
sorption features (e.g. atmospheric waterlines, etc.) we 
shifted the range by 200A and computed the median S /N 
for this range. The raw data were reduced using the 
HIRedux routine, then extracted, coadded and contin- 
uum fit with x_continuum. These routines are all part 
of the publicly available XIDL reduction package devel- 
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Table 1 

Journal of Observations 



QSO 


RA 


DEC 




Date Observed 


Exposure Time 




Av rcs b 


S/N c 




i 19000 V\\ 


( 19000 V\\ 






\ s ) 


(km s ) 




Q2359-02 


00 01 50.0 


-01 59 40.3 


2.800 


1997 Sep 29 


14400 


8 


17 










1997 Sep 30 


10800 


8 












1997 Oct 01 


12516 


8 




Q0000-2619 


00 03 22.9 


-26 03 16.8 


4.110 


1994 Sep 30 


10800 


8 


15 










1994 Oct 01 


3600 


8 




BR0019-15 


00 22 08.0 


-15 05 38.8 


4.530 


1996 Sep 20 


13500 


8 


18 










1996 Sep 21 


12600 


8 












1996 Sep 22 


14400 


8 




J0040-0915 


00 40 54.7 


-09 19 26.9 


4.976 


2011 Jan 16 


7200 


6 


7 










2011 Jan 24 


3600 


6 





Note. — (This tabic is available in its entirety in a machine-readable form in the online journal. A portion is 
shown here for guidance regarding its form and content.) 
a The emission redshift of the quasar. 

Av rcs is defined as the FWHM resolution of the spectrum. 
c The average signal to noise (S/N) ratio per 1.4 km s" 1 pixel. 



oped by J.X. Prochaska (jProchaska et al.ll2003bD . The 
reasons for choosing a single high resolution (i.e. R > 
40000) instrument are discussed in £12.31 

2.2. Methodology 

After the reduction process, we searched the reduced 
spectra for metal lines at the redshift of the DLA ob- 
tained from the Lyman alpha line in lower resolution 
data. The redshift of the DLA was then adjusted to 
coincide with the peak absorption feature in the metal 
lines, since metal lines are much narrower than the large 
width of the damped profile of the Lya line. This pro- 
vides a more accurate measurement of the DLA red- 
shift. We also calculated the column densities of met- 
als using the apparent o ptical depth method (AODM; 
iSavage fc Sembach 1991), which yields accurate column 
densities even if some of the lines are slightly saturated. 
From the column densities of the metal lines and the col- 
umn density of hydrogen, we can calculate the metallicity 
([M/H]) of the DLA, which is defined by: 

[M/H] = log 10 (N M /N H ) - log lo (NM/N H ) (1) 

The waveleng ths and o s cillato r strengths used in this pa- 
per are from IMortonl (12003D. wherea s the solar abun- 
dances are from lAsplund et all (|2009l) . 

Besides the metallicity we can also determine the cool- 
ing rate of neutral gas from these column densities. 
iPottasch et al.l ()1979[ ) showed that, i c , the cooling rate 
per H atom due to [C II] 158 /im emiss ion, which is 
the d ominant coolant in the neutral ISM ( Wrig hTet al.l 
is given by: 

N(C II*) 

4 = ^7 L A u ihv uh (2) 

where A u i and hv u i are the Einstein coefficient and en- 
ergy of the tra nsition from the ex cited to the ground state 
of C+ (see e.g. lWolfe et a l. 2003). We are able to measure 
the column density of C II*, because the C II*A1335.7 
fine structure line arises from the same state as the [C II] 
158 /im line, and falls within the spectral regime covered 
by optical telescopes for the redshifts examined in this 
paper. 

Note that an accurate H I column density, iVm, is re- 
quired to obtain accurate measurements of both metal- 



licity and £ c . We are generally unable to obtain the 
H I column density from the HIRES spectra, because 
accurate measurement of this quantity requires the spec- 
trum to be fluxed a nd HIRES spectra are difficult to flux 
(jSuzuki et al.l l2~003) . We therefore relied mainly on spec- 
tra obtained from the Sloan Digital Sky Survey (SDSS; 
Abazaiian eEaDHOQl) to measure Njj.i- However, for 
those 41 DLAs that were obser ved with the Echel ette 
and Imaging Spectrometer (ESI; iSheinis "eTaT][2002l) on 
the Keck II 10m telescope, we used the ESI spectra, be- 
cause the improved resolution provides more reliable and 
precise H I column densities. The column density of H I is 
determined from these spectra by simultaneously fitting 
a Voigt profile to the absorption profile of the Lya line of 
the DLA, and fitting the continuum of the fluxed spec- 
trum of the background qu asar (see e.g. iProchaska et al.l 
l2003at iRafeTski et al J 12012 ). 

The final two quantities discussed here are the two 
kinematic parameters, W1526 and Avgo. The reason for 
using two different kinematic parameters is that they 
describe different kinematical properties of the DLA. 
The Avg parameter describes the velocity width of the 
main neutral absorption complex by specifically ignor- 
ing weak outlying velocity components. On the other 
hand, the W1526 parameter is dominated by these out- 
liers, since this absorpt ion line is saturated in most cases 
(|Prochaska et al.ll2008D . Hence the W1526 parameter de- 
scribes the kinematics of the halo gas and/or weak satel- 
lite components. 

The rest equivalent width of the Si IIA1526 line, W1526, 
is defined as W1526 = W t s /(1 + 2), where W b s is the 
observed equivalent width of the Si IIA1526 transition. 
We choose the Si IIA1526 line because the line is mea- 
sured in the majority of spectra and its high oscillator 
strength causes the line to be saturated in most systems. 
It is important to stress that the W1526 statistic is almost 
independent of the column density of Si + ,N§[+ , if the ab- 
sorption line is saturated. Saturated lines are on the flat 
part of the curve of growth where the column density and 
equivalent width scale as, W1526OC -y/hi -^Si+ j anv reason- 
able change in would only marginally change Wi^q- 
There is therefore no a priori strong correlation expected 
between the equivalent width and any parameter which is 
derived from the column density of any metal lines, such 
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as metallicity, if the equivalent width is obtained from a 
saturated line. Consequently, we were mindful to only 
select those DLAs with saturated Si IIA1526 transition 
lines for the comparison of the W1526 parameter with any 
other parameter. However, the unsaturated lines were in- 
cluded in the distribution for W4526; otherwise we would 
bias this distribution towards higher equivalent widths, 
since low equivalent width systems are more likely to be 
unsaturated due to the mass-metallicity relationship. 

The velocity width, Avgo, of an absorption system is 
defined to be the width of the absorption profile in ve- 
locity space. To measure the veloci ty width of a DLA, 
we em ploy the same strategy as in Prochas ka fe Wolfd 
(|1997ft . with a few exceptions. We first select an un- 
saturated low-ion transition line. We require the line to 
be unsaturated, because a saturated line could overesti- 
mate the size of the velocity interval by including weak 
outlying velocity features that contain little of the to- 
tal neutral gas of the DLA. As such we require that the 
normalized flux, defined as F norm = I(v)/I c , where I c is 
the continuum intensity incident on the gas and I(y) the 
transmitted intensity, is greater than 0.1 over the entire 
absorption profile. A low-ion transition line is chosen 
because low-ions, such as Fe IIA1608 and Si IIA1808, are 
likely to trace the neutral gas, which creates the damped 
Lyma n a profile (jWolfe et alJ 119951 : IProchaska fe Wolfel 
1997). By contrast, higher ionization transition lines 
such as C IVAA1548, 1550 and Si IVAA1393, 1402 exhibit 
different velocit y structures (jWolfe fe Pr ochaskai 120001 : 
iFox et all 120071 ). In addition, we require that the line 
is unblended from any other absorption feature. 

After selecting the absorption profiles, we obtain an 
apparent optical depth profile; 

T(v)=\n[F Borm - 1 ]. (3) 

The resultant profile is smoothed to 8 km s _1 , the largest 
Av rcs of our sample, to prevent differences in resolution 
from affecting the Avgo values. After smoothing, we se- 
lect only those profiles for which the peak optical depth is 
detected at the 12c detection limit, so that components 
one-fourth this peak optical depth are detected with a 
3cr detection limit. This is important because these com- 
ponents could contain a significant fraction of the neu- 
tral gas content. T his criterion is less restric tive than 
the criterion used bv IProchaska fe Wo lfe (19973), because 
we believe that a 3er detection limit is enough to dis- 
cern small absorption features, since in almost all cases 
we have stronger absorption lines that clearly show the 
presence of these small absorption features above the 5a 
detection limit, and we are only after the velocity width 
of the line and not any of the other parameters described 
in IProchaska fe Wolfel (|1997j ). Finally, we find the width 
of the profile by stepping inward from both sides of the 
profile pixel-by-pixel, until we reach 5 percent of the to- 
tal integrated optical depth on each side. This width 
is the measured Avgo value. This last step prevents 
weak outlying absorption features from skewing the Avgo 
statistic. The complete set of smoothed optical depths 
as a function of relative velocity for all DLAs is shown 
in Figure [TJ For display purposes, we shift the data such 
that the left edge of the profile lines up with km s _1 . 
The smoothed l-er error array is shown as a (red) dashed 
dotted line. The separation between the (green) dashed 
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Figure 1. Optical depth profile as a function of relative velocity. 
The sample profile shown is of the DLA at redshift z a t, s = 2.9340 
towards QSO J1410+5111. The profile has been smoothed to 
8 km s" 1 , which is the largest Av res used in this sample. The 
dashed (green) lines mark the Avgo interval, which is 243 km s _1 , 
and the dashed-dotted (red) line is the smoothed 1-ct uncertainty 
on the data. For the whole sample, the Avgo values range between 
17 km s" 1 and 484 km s -1 . 

lines marks the velocity width. We assume an error of 
10 km s -1 on these meas urements, similar to that used 
bv IProchaska et all (|2008f ). 

2.3. Systematic Errors and Biases due to Sample 
Selection Criteria 

There are two main concerns in selecting an appro- 
priate sample for exploring correlations: systematic er- 
rors and biases caused by sample selection criteria. One 
source of systematic error could be the use of two dif- 
ferent resolution instruments, which would in particular 
affect the velocity width measurement since it is mea- 
sured by moving in pixel-by-pixel in the spectra. Lower 
resolution spectra, like the spectra taken with ESI which 
has a resolution of 44 km s _1 , would overestim a te the 
velocity width. This motivated IProchaska et al.l (2008) 
to reduce their ESI Avgo measurements by 20 km s _1 . 
To circumvent the systematic errors caused by multiple 
instruments, we decided to restrict ourselves to veloc- 
ity width measurements and metal column densities ob- 
tained from HIRES alone. 

The other main concern is the effect of sample selection 
biases on the data. To assess this effect, we recorded the 
reason why each DLA was observed by our group with 
HIRES. The reason generally falls in one of the follow- 
ing four categories: (1) H I-selected, the DLA was se- 
lected based on just the column density of H I; (2) high 
redshift, the DLA was selected because of its high red- 
shift; (3) metal content, the DLA exhibited very strong 
metal lines in its spectrum obtained from SDSS; and fi- 
nally, (4) serendipitous, the DLA happened to fall in the 
same line-of-sight of the primary target. In addition, 
the majority of targets were observed because the flux 
of the background quasar exceeded a certain threshold; 
although the second and third reasons often drove our 
group to observe fainter targets. Each of these selection 
criteria needs to be examined for any bias that the corre- 
sponding subsample might contain. First, if the distance 
between the DLA and quasar is sufficient, then the DLA 
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Table 2 

HIRES DLA sample 



QSO 2 abs logiVui Metallicity log £ c W1526 Avgo a Transition 13 Selection References 

(cm -2 ) (ergs s -1 H — 1 ) (A) (km s _1 ) Criterion 



Q1104-18 


1.6613 


20.80 ± 0.10 


-0.99 ± 0.10 


-26.87±0.11 




50 


Si 


IIA1808 


H 1 - selected 


6 


Q1331+17 


1.7763 


21.14 ± 0.08 


-1.34 ± 0.08 


-26.60±0.09 


0.499±0.001 


72 


Si 


IIA1808 


H 1 - selected 


7, 13, 18, 24 


Q0841+12 


1.8640 


21.00 ± 0.10 


-1.46 ± 0.10 






30 


Si 


IIA1808 


Serendipitous 


20 


Q2230+02 


1.8643 


20.85 ± 0.08 


-0.71 ± 0.09 




1.500±0.018 


172 


Si 


IIA1808 


H 1 - selected 


7. 13 


Q1210+17 


1.8917 


20.60 ± 0.10 


-0.79 ± 0.10 




0.419±0.003 


38 


Si 


IIA1808 


H 1 - selected 


13, 19 



References. — (1) Wolfe ct al. (1994); (2) Lu ot al. (1996);(3) Prochaska & Wolfe (1996);(4) Prochaska & Wolfe (1997); (5) Lu ct al. (1998); (6) 
Lopez ct al. (1999); (7) Prochaska & Wolfe (1999); (8) Molaro ct al. (2000); (9) Petitican ct al. (2000); (10) Prochaska & Wolfe (2000); (11) Ellison et al 
12001al); 7l 2) IMolaro et all ( 120011); (13) IProchaska et al.1 <2001l); ( 1 4) | Ledoux et al.l d2002 | ); ( 15) ILevshakov et al.l I2002 | ); (16) IProchaska et al 



(20) ILedouxetal 



(25) IVladilo et al 



(2006 
2011 



(21) 
(26) 



(17) Prochaska ct al. (2003a); (18) Dcssaugcs-Zavadsky ct al. (2004); (19) Dcssaugcs-Zavadsky ct al. (2006) 
IDessauges-Zayadskv et al.l <2007l) ; (22) IProchaska et al.l 120071) ; (23) IWolfe et all 120081) ; (24) IJorgenson et all ( 120101) 
IRafelski et alj (120121 ) 

Note. — (This tabic is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its 
form and content.) 

a The uncertainty on the Avgo measurement is taken to be 10 km s — 1 for each measurement, 
k Transition used to measure the Avgo parameter. 



properties are unrelated to the properties of the back- 
ground quasar. This indicates that selecting DLAs by 
the flux of the background quasar should not bias the 
DLA sample. Similarly for the serendipitous subsam- 
ple, if the separation between the primary target and 
the serendipitous DLA is large enough, the serendipitous 
DLA should be unrelated to the primary target. We take 
a velocity separation of 3000 km s _1 to be sufficient for a 
DLA to be unaffec ted by another DLA or quasar (see e.g. 
Ellison et alll2010l ). and we have made sure that no DLA 
in our sample violates this criterion. Secondly, the DLAs 
that were observed because of t heir metal column density 
(metal-strong DLAs, MSDLAs; [Herbert Fort et al.ll2006[) 
are not included in our sample. This subsample is biased 
in metallicity and W1526, because metal-selected DLAs 
have on average high er metallicities and correspondingly 
higher W1526 values ([Kaplan et al.l 120 lOh . 

This leaves two sample selection criteria; the redshift 
and the H I column density of the DLA; both of these 
sample selection criteria cause a bias in the sample. The 
bias in the redshift is twofold. First, bright quasars are 
more abundant in the redshift range 2-3, and therefore 
any magnitude-limited selection of qua sars will produce 
a ma j ority of DLAs within thi s range ([Prochaska et al.l 
[20051 IProchaska fc Wolfe! l2l)09h . Second, for a subset of 
the data we selected the quasars based on their high red- 
shift, which biases that sample towards DLAs with high 
redshift. The H I column density selection criteria causes 
our sample to contain proportionally too few DLAs below 
a column density of 10 20 5 cm~ 2 compared t o the sample 
of DLAs from SPSS. IProchaska et all (|2007t ) ascribe this 
to our preference in selecting higher column density sys- 
tems to ensure the absorber satisfies the DLA H I column 
density criterion. 

To gauge the extent of the bias on the distribution 
of the remaining four parameters; we want to compare 
our biased sample to one that is 'free' of the bias caused 
by the redshift and H I column density selection crite- 
ria. To accomplish this, we randomly divide our sam- 
ple of 100 DLAs in half. One half is left untouched, 
whereas we randomly pick 50 DLAs (with repeats) from 
the second half such that the redshift and H I col- 
umn density of the second half reproduces the H I col- 



umn density frequency distribution observed in DLAs, 
(hAN.X): IProchaska eFall [20051 [ Prochaska fc Wolfd 
2009; iNoterdaeme et al.l 120091 12012ft . This is indicative 
of a sample that is 'free' of biases in redshift and H I 
column density. To check for a bias in metallicity, Avgo , 
W1526! and t c of the untouched half, we look at the dis- 
tribution that these parameters assume for both halves 
of the data set. We compare the distributions of each 
half, using a two-sided Kolmogorov-Smirnov test (KS) 
test and a Mann- Whitney U-test. The KS test provides 
a probability that both populations are drawn from the 
same parent population, whereas the U-test provides a 
probability that the two median values are significantly 
different. This procedure of randomly dividing the sam- 
ple, reselecting one half of the sample and comparing the 
distribution is repeated 1000 times, to get a median value 
for the results of both tests. 

The resultant median values of the KS test for all pa- 
rameters are greater than 0.05 indicating that we cannot 
rule out the null hypothesis, which is that the two sub- 
samples were drawn from the same parent population, at 
a 95 % confidence level (c.l.). Similarly, the U-test yields 
mean values greater than 0.05 for all parameters indicat- 
ing that we also cannot reject this null hypothesis, which 
is that the medians of the two subsamples are the same. 
Hence, we see no evidence that the redshift and H I col- 
umn density selection criteria have significantly affected 
the distribution of the remaining parameters. It is im- 
portant to note that this does not mean that the redshift 
and H I column density show no correlation with the re- 
maining parameters. It means that the effect of the two 
sample selection criteria does not significantly affect the 
distributions of the remaining four parameters. Together 
with the serendipitous sample, the H I-selected and red- 
shift selected subsample comprises the complete sample 
used in this paper. In total, it consists out of 100 DLAs 
for which we have accurate measurements of metallicity, 
redshift, H I column density, and Avgo; the full sample 
is shown in Tabled 

3. DISTRIBUTIONS AND CORRELATIONS 

In this section, we compare the distribution of metallic- 
ity, £ c , and IU1526 parameters to previous studies specif- 
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Table 3 

Fitting Functions to the Distributions 



x a\ fj.2 o"2 r 



log Av 90 1.83 0.39 

log W1526 -0.44 0.42 

[M/H] -1.46 0.55 — 

log £ c -27 .4 0.10 -26.7 0.28 0.55 



Function Used 

( x M ) 2 

(a) f(xin u o*) = - r i—e ^~ 



(b) f(x; m,a?) = (j^)f(x; (ti,af) + (j^)f(x; ^^\) 



that the distribution is reasonably well described by a 
Gaussian with a mean metallicity of -1.57 and a disper- 
sion of 0.57. Our smaller data set with a mean metal- 
licity of -1.46 and dispersion of 0.55 is consistent with 
their sample; a one-sided KS-test shows that the smaller 
second peak at -0.61 is not statistically significant. Sec- 
ondly, the £ c distribution shown in Figure [2b, was first 
described by I Wolfe et"ail (|2008D ; they found that their 
sample of 37 positive detections of £ c showed a bimodal 
distribution. Our sample contains 32 positive detections, 
and shows a similar bimodality with peaks at 10~ 27 ' 4 ergs 
s _1 H _1 and 10 -267 ergs s _1 H _1 respectively. This al- 
lows the sample to be divided into two subsamples: a 
high-cool subsample, where £ c > 10~ 27 ergs s _1 H _1 , 
and a low-cool subsample for which £ c < 10~ 27 ergs 
s -i h~i Final l y, the W1526 distribution, discussed in 
IProchaska et a l. (2008), is shown in Figure [3]:. While 
most DLAs exhibit saturated Si IIA1526 lines, this distri- 
bution includes DLAs with unsaturated Si IIA1526 lines 
to prevent a bias against low equivalent widths (see Sec- 
tion [2~3]). The distribution is ful l y con sistent with the 
sample found in IProchaska et aTl (20081 ). The fact that 
there exists good agreement between the sample consid- 
ered in this paper and the larger samples used to study 
these three parameters in detail, shows that the current 
sample is a representative sample of DLAs for each of 
these three parameters. This is crucial for a proper ex- 
ploration of correlations between the parameters. 

The last parameter we consider is the Avgo statistic. 
The distributi on is shown in Figure J2}i. Unlike pre- 
vious papers (jProchaska fc Wolfdll^97nprochaska et all 
2008), we do not use linear bins, but divide the bins log- 
arithmically. This choice was made because for the cur- 
rent sample size, the velocity width distribution is well- 
approximated by a log-normal distribution. The Avgo 
distribution has two important characteristics. The first 
important characteristic of this distribution is the sharp 
decline with decreasing Avgo of DLAs with AV90 below 
17 km s _1 (log Avgo= 1.23). This decline is not likely to 
be an observational limit, because even our lowest reso- 
lution data has a Av rcs of 8 km s^ 1 which corresponds 
to a Avgo of 11 km s _1 for a single Gaussian component. 
Although some objects have been fou nd with these small 
Av 9 q values, they are extremely rare (iCooke et aTl 120111 : 
iPettini et all 120081; iPenprase etall 120101) . The rarity of 
these objects suggests that the majority of DLAs have at 
least one component with an internal velocity dispersion, 
a, greater than 7 km s _1 or have multiple components 
which exhibit velocities of ~ 10 km s _1 relative to each 
other. 

The second characteristic of the distribution of ve- 
locity widths is the large number of high Avgo DLAs. 
When plotted in a linear histogra m, these DLAs form 
a 'tail ' in the distribution, as s een in IProchaska fc Wolfe! 
(|1997ft : lProchaska et al.l(|2008l) . Most of these large Av 90 
DLAs have Avgo values below 200 km s -1 , but there are 
a few DLAs which show a velocity widt h significantly 
above this value (e.g. lPetitiean ct al. 2002J). These DLAs 
with very large velocity width are likely due to certain 
quasar sightlines encou ntering multiple galaxies or winds 
(jProchaska et al.l [200 8) . As a result of the large num- 
ber of high Avgo DLAs, the resultant median Av go is 
significantly higher than expected from numerical simu- 



Note. — (a) is a normalized Gaussian distribution function with 
mean, fj,\ and variance erf. (b) is the sum of two Gaussian func- 
tions. The factors in front of the Gaussian terms are required such 
that the total integrated area under the function is equal to unity. 
The r-parameter is the ratio of the relative sizes of the two Gaussian 
distributions. 

ically aimed at exploring these parameters, to show that 
our sample of 100 DLAs is a representative sample of 
each of the parameters. We explore our sample of 100 
Avgo measurements of DLAs, which to our knowledge 
is the most accurate distribution to date of this quan- 
tity. We also look at any potential correlations between 
the parameters, except for those that have been explored 
in detail before, such as the kinematics-metallicity cor- 
relation (iWolfe fc Prochaskal 119981: iLedoux et all 120061 : 
IProchaska et al. 120081) and the metallicity-redshi f t cor- 
relation (IProchaska et~aT1 12003a! : iKulkarni et~aT1 120051 : 
iRafelski et al.ll2012l ). 

3.1. Distributions 

The distributions of metallicity, £ c , and W1526 of the 
sample are shown in Figure [2ja-c). These distributions 
have been discussed in detail in previous papers and 
here we give a short summary of the characteristics of 
each distribution. First, we scale each distribution such 
that its integrated area is normalized to unity. Conse- 
quently, we can interpret the y-axis of parameter, p, as 
the probability of finding a DLA within the parameter 
range (p,p+dp). We then fit analytic functions to each of 
the distributions described above; these are shown as the 
solid (red) line in Figure [2] To prevent binning from af- 
fecting the fit, we do not fit the actual distributions, but 
instead fit the cumulative distribution function (shown 
as insets in Figure [2]) . The procedure used for the fitting 
process is as follows. We first use a chi-squared fitting 
technique to fit the cumulative distribution function of a 
single Gaussian distribution function to the cumulative 
distribution of the data; this fit is then compared using 
a one-sided KS-test. If the fit is unacceptable, we try for 
the fitting function a sum of two Gaussian distributions. 
The resultant equations to the fit are shown in Table [H 
As expected, the £ c sta tistic is best desc ribed by a bi- 
modal fitting function (W olfe et al.l l2008h . whereas the 
other parameters can be described by a single Gaussian 
distribution. 

Figure^ shows the metallicity distributio n, which was 
recently published by Rafclski et al. (2012j). They show 
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Figure 2. Distributions of metallicity, 158 fim cooling rate, equivalent width and velocity width of the sample of DLAs. The solid (red) 
line is a chi-squared best fit line to the cumulative distributions (shown as insets in the figure); the equations of these best fit parameters 
are given in Table [3] The distributions of metallicity, Avgo and W1526 can all be described by a single Gaussian distribution. The i c 
parameter on the other hand, is clearly bimodal and is best described by a bimodal distribution. The distributions of metallicity, £ c and 
W1526 are in good agreement with previous studies, indicating that the sample is a representative subsample of each of these parameters. 
The Avgo distribution shows a large number of DLAs with Avgo > 100 km s _1 , which raises the median Avgo value to 70^jS km s — 1 . 



lations (iProchaska &: Wfill997tlRazoumov et al.ll2006t 
iPontzen et al.l 20081 ). The median value for our sample 
is 701^3 km s _1 . The uncertainty on the median is the 
combined uncertainty of three different factors. First, 
due to the finite sampling of the Avgo distribution, our 
sample of 100 DLAs might not evenly be divided by the 
true median of the sample. To estimate the uncertainty 
due to this factor, we note that this is a simple selection 
problem which can be described by the binomial distri- 
bution with p = 0.5. The binomial distribution has a 1-a 
uncertainty of 1/2 • y/n, where n is the number of data 
points. This means that for a sample of 100 DLAs, there 
is a probability of 0.68 that the true median of AV90 is 
contained in the range between DLAs 45 and 55 when 
the DLAs are arranged in order of increasing Avgo. To 
this uncertainty, we need to add in quadrature the un- 
certainty arising from the limited sampling of the data. 
This uncertainty we take to be the average spacing of the 
y/n points centered around the median. Finally, we add 
in quadrature the mean observational uncertainty to get 
the uncertainty on the median. 



3.2. Correlations and Dependencies 

In this section we discuss some of the two-parameter 
correlations and dependencies that exist between the six 
parameters discussed in this paper. Table 0] lists all 15 
combinations of two parameters. For all 15 combinations 
of parameters, we are not only interested in any potential 
correlation, but we also want to know if the distribution 
of one parameter is dependent on the other. Therefore we 
apply a variety of different tests described below which 
test for both correlations and such dependencies. 

To test for the existence of any potential correlation, 
we fi nd a linear fit to the data using the MPFITEXY rou- 
tine dWilliams et al.ll2010D. wh ich depends on the MPFIT 
package by IMarkwardtl (|2009|) . This routine takes into 
account the uncertainties in both parameters when cal- 
culating the slope and y-intercept of the best fit line. To 
calculate the uncertainty of the slope and y-intercept, we 
use a bootstrap method. The bootstrap method works 
by randomly selecting 100 DLAs from the original sam- 
ple of 100 DLAs, but allowing for repeats. The resultant 
sample is then fitted using the same fitting routine as 
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Table 4 

Table of Potential Correlations 



Dependency Linear regression line a Kendall b KS-test b F-test b U-test b T-test b 
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a These are the best fit parameters for the linear regression line of the form: y — ax + b 
b Values are probabilities that the null hypothesis of each test can be rejected 
c Not enough data points are available to get an accurate measurement 



the original sample, and the slope and y-intercept are 
recorded. This sampling and fitting is repeated 1000 
times, creating a distribution of slopes and y-intercepts. 
The 1-cr uncertainty on the slope and y-intcrccpt is in- 
ferred from a Gaussian fit to these distributions. A sec- 
ond test for the existence of a correlation between the 
parameters is provided by the application of a Kendall 
Tau test. The resultant two-sided significance of its de- 
viation from zero is shown in Table |4j here a small value 
indicates a probable correlation. 

To test if the distribution of one parameter is depen- 
dent on a second, we apportion our complete sample into 
three equally sized subsamples based on the value of the 
second parameter. We then compare the distribution of 
the first parameter for the subsamples with the small- 
est and largest values of the second parameter using four 
different tests. The first test we apply is a two-sided 
Kolmogorov-Smirnov (KS) test. This test will provide 
a probability that the two subsamples are drawn from 
the same parent population. A value smaller than 0.05 
indicates that we can assume that the two populations 
were drawn from a different parent population at the 
95% confidence level (c.l.). To test if the variance of the 
two subsamples is significantly different, we use the F- 
statistic. Again a value of 0.05 or smaller indicates a 
significantly different variance in the two subsamples at 
the 95% c.l. Finally, we used a Student's T-test and the 
Mann- Whitney U-test to compare the mean and median 
of the two subsamples, where again a smaller than 0.05 
probability would indicate that the two subsamples have 
significantly different means or medians. The results of 
these tests for all 15 combinations are shown in Table 01 

In the next few subsections we will discuss some possi- 
ble correlations that have not been discussed in previous 
papers; in particular the dependence of the H I column 
density to the other parameters. We also explore the pos- 
sible correlation between redshift and kinematics, which 
is important for exploring the interplay between redshift, 
metallicity and mass of DLAs. 



3.2.1. [M/H] vs Nhi 

Figure [3] shows a plot of metallicity versus H I col- 
umn density for the objects in our sample. The linear 
line was calculated using the method described in the 
previous section, and the shaded gray area marks the 1- 
<7 uncertainty on this line. There are two features of 
this plot we would like to discuss. First, the sample 
shows a lack of high H I column density - high metallic- 
ity systems. IBoisse et~a l. (1998) attributed this lack to 
a dust bias; presumably these systems would have large 
enough dust fractions to block out all of the light from the 
QSO and therefore escape detection. This assertion has 
been question e d by many pape r s incl uding Elliso n et al. | 
(|2001bL l2005h: lAkerman et al.[ (|2005[ k Uorgenson et al.l 
P006f k iFrank 



Perouxl ( 201C ) who found that radio- 
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Figure 3. Metallicity versus H I column density. The square 
points mark our sample, and the solid line is a linear fit to the data, 
where the gray area marks the 1-cr error on this line. The sample 
shows a lack of high H I column density - high metallicity systems. 
T he lack of the s e syst ems has been noted in previous papers such 
as IBoisse et al.l 119981V However, the sample also shows a lack of 
high H I column density - low metallicity systems. Using the F- 
statistic, we show that the lack of these systems is not likely (i.e. 
less than 5% probability) due to the small number of high H I 
column density systems. 
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Figure 4. Remaining four parameters plotted against the H I column density. The figure markings are the same as Figure(3] Both redshift 
and Avgo are consistent with no correlation, whereas W1526 shows a potential correlation with the H I column density, although this is 
not statistically significant for the sample size in this paper. The l c distribution tentatively shows that the bimodality is more evident at 
higher column densities. The dotted line indicates the line below which we cannot reasonably positively detect a cooling rate, because of 
the required S/N needed for a positive detection. 



selected DLAs do no t differ signifi c antly from those se- 
lected optically, and IKaplan et al.l (|2010f) who detected 
certain meta l -stron g DLAs above the threshold found in 
iBoisse et al.l (1 19981). A second in terpretation for this lack 
was put f orth by iSchavd (120011). which was further ex- 
plored by Krumho lz et al.l (|2009f ). who showed that the 
absence of high H I column density systems could be due 
to a transition from the atomic to molecular state of the 
atoms in the cold phase of a two-phase medium. How- 
ever, both of these explanations only describe the lack of 
high column density - high metallicity systems, and are 
unable to explain the second feature, which is the lack of 
high H I column density - low metallicity systems. In 
fact, both features seem to be quite symmetric, in that 
the scatter plot exhibits a reflection symmetry about the 
line [M/H] = —1.43. To test if these features are statis- 
tically significant, we use the F-statistic. Table [4] shows 
that the F-statistic for the [M/H] vs TVm correlation is 
0.013. This indicates that the null hypothesis, which is 
that the smallest H I column density DLAs (i.e. DLAs 
with A^hi < 20.5) and the largest H I column density 
DLAs (i.e. DLAs with Am > 20.85) have the same vari- 
ance, can be ruled out at a 95% confidence level. Hence, 
both features are not likely due to small number statis- 
tics; possible explanations for the existence of these two 



features and the symmetry between them are given in 
Section ISTTI 

3.2.2. Nhi dependencies of the remaining parameters 

Figure 2] shows the remaining four parameters plotted 
vs H I column density. Both redshift and Avgo show 
no statistically significant correlation with JVm, and the 
distribution of these two parameters is also not depen- 
dent on Nhj. For VF1526 there is a greater likelihood that 
we see evidence for an increase in equivalent width with 
H I column density as indicated by the Kendall tau pa- 
rameter. However, this trend is not significant. This can 
be seen in the slope of the linear regression line which is 
consistent with no correlation at the 2-cr level. Finally, 
for the £ c parameter, the high Am subsample has a sig- 
nificantly smaller mean and median than the low Am 
subsample. This is due to two effects. First, we can ob- 
serve a similar trend as in the above case for metallicity 
that higher H I column density systems tend to have less 
variance compared to lower H I column density systems. 
Since the £ c distribution is bimodal, the higher H I col- 
umn density systems have values close to either of the two 
mean values, and therefore show a more distinct bimodal- 
ity than the lower H I column density systems. Secondly, 
we are unable to measure low H I column density, low- 
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Figure 5. Redshift evolution of the two different kinematic parameters. The AV90 statistic shows no evidence for any redshift evolution. 
The linear fit is clearly consistent with a flat line. The W1526 statistic shows potential evidence for redshift evolution as can be seen by the 
solid linear regression line. However, the sample size is too small to confirm this correlation. 



cool systems, because we only use positive detections of 
the C II*A1335.7 line for calculating £ c , and due to the 
limited signal to noise ratio of the spectra we are unable 
to positively detect the very low CII* column density 
systems. Both these effects contribute to the differing £ c 
distribution for the low Nn I and high TVhi subsamples. 

We would like to point out that this second effect also 
biases the relative size s of the low-cool and high-cool sub- 
samples. Wolf e et al.l (|2008l ) report relative sizes of the 
two subsamples to be 38% for the low-cool subsample 
and 62% for the high-cool sample using only positive de- 
tections. To gauge the extent of this bias, we estimate 
that for a spectrum with a signal to noise ratio of 30, the 
minimum value we can measure for the column density 
of CII* is approximately 10 125 cm~ 2 . Using, equation [5] 
we can compute the corresponding £ c value as a function 
of H I column density. This boundary, below which we 
cannot make any positive detections, is shown in Figure 
Hi by a dotted line . This shows that the result quoted in 
Wolfe et al.1 ((2008) underestimates the number of DLAs 
in the low-cool sample. Comparing the number of DLAs 
in the two subsamples for DLAs with H I column density 
above 10 20 ' 7 cm -2 , for which this bias is negligible, shows 
that there are in fact more DLAs in the low-cool subsam- 
ple than the high-cool subsample. However, we caution 
that not enough data points are available to make a pre- 
cise estimate of the sizes of the two subsamples. 

3.2.3. Avgo and W152G vs redshift 

The last two possible correlations that we consider are 
the two kinematic properties, Avgo and W1526, as a func- 
tion of redshift. The kinematical properties are believed 
to be strongly linked to the mass of the dark matter 
halo. T his is a ubiquitous feature of almost all DLA mod - 
els (e.g. IHaehnelt et al.l[l998t IProchaska fc WolfdH997l ). 
Since current ACDM models predict hierarchical galaxy 
formation, which in turn predict an increase in the mass 
of dark matter halos over time, we may expect to see a 
corresponding increase in the Avgo and 1^1526 statistic 
of the DLAs as well. Again we use the tests described in 
Section [3~2l to see if there is any linear trend in the data, 
and/or if the distribution of Avgo or W1526 is evolving 



with redshift. 

Figure [5^ shows Avgo as a function of redshift. As in 
Figure 02 the solid black line indicates the best fit lin- 
ear line to the data, and the gray area marks the 1-a 
uncertainty on the parameters. The resultant best fit 
is: logio(Av 90 ) = (0.00±0.04)-z + (1.84±0.13). This lin- 
ear regression line is clearly consistent with no redshift 
evolution. Furthermore, all other tests show that the 
distribution of Avgo is also consistent with no evolution. 

This is in c ontra st with a previous result found by 
iLedoux et al.1 (|2006l ). who found that the mean Avgo de- 
creased with increasing redshift when they apportioned 
their sample into two subsamples based on redshift. 
When we apportion our sample into two subsamples 
based on redshift, we find two main differences between 
their results and ours. First of all, our sample covers 
a larger redshift range; in particular our median values 
for each redshift subsample are 2.400 and 3.722 whereas 
their median redshifts are 2.087 and 2.796; any redshift 
evolution should therefore be more noticeable in our sub- 
sample. Second of all, their sample was smaller and no 
errors were reported in the median value of the velocity 
widths. Indeed, if we calculate the errors on the median 
for their two subsamples, we get a median of 69l 2 3km s _1 
and 89 t 2 g for the high redshift and low redshift subsam- 
ple respectively. This shows that the two medians are 
within 1-a of each other, indicating that their sample is 
also consistent with no redshift evolution. 

The second kinematic parameter we consider is W1526, 
shown in Figure [5]d. The data shows a lack of high- 
rcdshift large-equivalent width systems, and the Kendall 
Tau correlation test shows a small value of 0.042 indicat- 
ing that a potential correlation could exist. However, a 
linear fit to the data gives logio(W r i526) = (— 0.10±0.05) • 
z — (0.08±0.15), which indicates that the correlation is 
seen at only the 2-cr level. Moreover, when we compare 
the median and mean of the highest redshift DLAs with 
those of the lowest redshift DLAs, we do not see a signif- 
icant difference in their value. Indeed the KS test shows 
a likelihood of 56% that the two samples are drawn from 
the same parent population. We therefore conclude that 
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Figure 6. Panel (a) of this figure illustrates that when redshift, Avgo and metallicity of a DLA are plotted on a three dimensional scatter 
plot, the points trace a plane in this parameter space, although with scatter. This plane is also seen when the AV90 parameter is replaced 
by Wi$26 (panel (b)). The plane is marked by thick, dark gray lines. The data points are connected to this plane by solid black (gray) 
lines if they fall above (below) the plane. The grayscale (color scale) of the data points is correlated to the metallicity of the DLA. 
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Figure 7. This figure illustrates the advantage of using the fundamental plane equation in describing the two correlations shown in 
the top two figures. Panel (a) shows the metallicity-Avgo correlation where we have apportioned our sample into two equal subsamples 
based on redshift. The solid (green) triangles are the higher redshift DLAs, whereas the solid (purple) squares are lower redshift DLAs. 
This panel clearly shows that a large part of the scatter is due to the correlation between redshift and metallicity. Similarly, the scatter 
in the metallicity-redshift correlation (panel (b)) is due to the existence of a mass-metallicity correlation at each redshift. By applying 
the fundamental plane equation we can reduce this scatter as is shown in panels (c) and (d). The scatter is still bigger than the mean 
observational uncertainty, which is shown by the errorbar. 



with the current sample size, the correlation between 
Wi52e and redshift cannot be determined at > 3-<r sig- 
nificance level. 

4. THE FUNDAMENTAL PLANE(S) OF DLAS 

As was discussed in the introduction, the aim of this 
paper is to examine the interplay between multiple pa- 
rameters of DLAs using the fundamental plane descrip- 
tion. For this description to be useful, we need to have 
three parameters where two parameters are not strongly 
correlated, but each parameter shows a strong correlation 
to the third parameter. In Table 2] we list all of the corre- 
lations that exist between the six parameters discussed in 
this paper. The only parameters that clearly satisfy this 
criterion are metallicity, redshift and Avgo, and metal- 
licity, redshift and VFi.526. When we plot the first three 
parameters for each DLA on a three-dimensional scatter 
plot, we see that the points indeed fall close to a plane 
inside this space, although with scatter (Figure [5^)- Fig- 
ures [T^i and[7jD are two projections of this plane along the 
redshift and Avgo axis respectively, where the third pa- 
rameter is apportioned into two eq ually sized subgroups . 
Figure [7^ is similar to Figure 2 of ILedoux et~aTl (12006D 
and F igure [TJd is similar to Figure 11 of iRafelski et all 
(2012). Figure [7^, shows that the scatter in the Avgo- 
metallicity correlation is in part due to redshift evolution 
of metallicity and similarly the scatter in the redshift- 



metallicity correlation is in part due to differences in the 
kinematics of the DLAs at each redshift. 

To accurately describe the shape of the distribution in 
this three-dimensional space, we need to fit a plane equa- 
tion to the scatter plot. To do this, we us e the 'direct fit' 
least s quare fitting technique described in Bcrnard i et al.l 
(2003), which fits a fundamental plane equation of the 
form [M/H] = a-logAvgo + b-z + c to the data set. To 
account for measurement errors, we weigh each individ- 
ual point by the inverse of the estimated uncerta i nty us - 
ing the methods described in ISheth fc Bernardil (2012). 
This produces the following covariance matrix, R which 
is normalized by the RMS value of each quantity: 

Av 90 z [M/H] 

Av 90 / 1.000 0.121 0.361 \ 

z 0.121 1.000 -0.396 

[M/H] V 0.361 -0.396 1.000 / 

From this matrix we can calculate th e coefficients of 
the pl ane using the equations described in lBernardi et all 
( 2003) . The resultant plane is described by the following 
equation: 

[M/H] = (-1.9±0.5) + (0.74±0.21)-logAv 90 

-(0.32±0.06)-z ^ 

The uncertainty of the parameters are estimated using 
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Figure 8. This figure illustrates the fundamental plane using the W1526 statistic instead of Avgo. The panels are the same as Figure [71 
except the kinematic parameter has been switched. Note that the scatter in the equivalent width-metallicity relationship is not significantly 
reduced by applying the fundamental plane equation. This is due to the fact that the plane is only slightly tilted in the redshift direction. 
This, in part, explains the observational result t hat the equivalent width-metallicity relationship exhibits the smallest scatter of any 
parameter with metallicity (jProchaska ct al. 20QH). 



the bootstrap method discussed in Section [3~2l It is im- 
portant to note that because Avgo and redshift are not 
correlated (Figure [S£i), the uncertainties in the coeffi- 
cients of the fundamental plane are independent of each 
other. Figure [7J: and [TJi are plots in which we take into 
account the third parameter using the plane equation de- 
scribed above (i.e. this is like rotating the plane until it 
is edge-on). To quantify the reduction in scatter, we cal- 
culate the RMS scatter around a linear fit to the scatter 
plot before and after we apply the plane. The scatter 
in the correlations before we apply the plane equation 
is 0.45 dex (0.50 dex) for the metallicity-Avgo (redshift- 
metallicity) relationship. After we apply the plane the 
scatter is reduced by approximately 20 % to 0.37 dex 
(0.38 dex) for each of the correlations. The reduction in 
scatter is clearly visible in Figures [TJ: and Hi. However, 
the scatter is still significantly bigger than the observa- 
tional uncertainty on each measurement, which is on av- 
erage 0.12 dex. We will comment more on the reduction 
in scatter in Section [521 

The second fundamental plane of interest is simi- 
lar to the first, except the Av go statistic is replaced 
by the other kin e matic parameter, W1526 (Figure |BJd). 
iProchaska et al.l ((2008) showed that the metallicity- 
W1526 correlation exhibits smaller scatter than the 
metallicity-Avgo relationship, and a fundamental plane 
equation between metallicity, redshift and W1526 may re- 
duce this scatter further. As in Sections 13.2.21 and 13.2.31 



we are using just the subset of W1526 measurements with 
saturated W1526 values, although the inclusion of these 
systems does not significantly affect the shape or tilt of 
the fundamental plane. The equation of the fundamental 
plane between redshift, metallicity and W1526 is: 

[M/H] = (-0.47±0.14) + (l.l±0.3)4ogWi52 6 

-(0.16±0.06)-z ' 

Figure [8^, shows metallicity versus W1526', unlike Fig- 
ure [7^ where the Avgrj-metallicity trend moves down- 
ward with increasing redshift (i.e. at higher redshift, 
a given Avgo corresponds to a lower metallicity), the 
W / i526 _m etallicity correlation does not evolve with red- 
shift. The lower redshifts are higher on the correlation 
trend line because of the correlation between redshift and 
metallicity, but very little of the scatter is due to redshift, 
and therefore a fundamental plane description does not 
reduce the scatter significantly (Figure [5J;) To be spe- 
cific, the scatter is reduced by only 0.02 dex. Visually 
this means that we are already looking edge-on to the 
plane in Figure [SJi. However, this second fundamental 
plane is able to reduce the scatter significantly (0.15 dex) 
in the redshift-metallicity relation (Figure [8jl) as was the 
case with the previous fundamental plane. The measured 
scatter around the two correlations after applying this 
fundamental plane is 0.29 dex for the Wi526-metallicity 
correlation and 0.30 dex for the redshift-metallicity cor- 
relation. The scatter around this plane is significantly 
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smaller than the scatter around the fundamental plane 
found previously, which we discuss further in Section f5.2l 

5. DISCUSSION 

The primary aim of this paper has been to explore the 
fundamental relations that exist between the parameters 
that describe DLAs. In this section we discuss the results 
of this study. 

5.1. Implications of two parameter dependencies 

Table |4] lists all of the possible two parameter depen- 
dencies for DLAs that are explored in this paper. The 
first notable dependency is that between the H I column 
density and the m etallicity of the DL A (Section I3.2.1|) . 
Previous work (e.g. Boisse ct al. 1998) observed a lack of 
high H I column density, high metallicity systems. Our 
sample, which was selected in a less biased way, shows 
that this lack of high H I column density, high metallic- 
ity systems is just half of the picture, and there is also a 
lack of high H I column density low metallicity systems. 
In Section GOT] we have given two possible explanations 
for the lack of high H I column density, high metallicity 
systems, but we need an explanation for the lack of high 
H I column density, low metallicity systems. 

One possibility is that this lower 'envelope' is not re- 
lated to the upper 'envelope'. The existence of just a 
lower envelope can be explained by the following two 
explanations: (1) higher H I column density DLAs are 
due to sightlines probing the inner part of g alaxies, be- 
cause of the metallicity gradients in DLAs (jChen et al.l 
2005), these sightlines should increase in metallicity. (2) 
Low metallicity DLAs are due to sightlines crossing cold 
flows, which because of their low density and high ion- 
ization level are unable to produce sightlines with a high 
H I column density. The latter is see n in si mulations 
modeling cold flows by Fumagal li et al.l (|201lD who find 
that only at higher redshifts cold flows are able to host 
low H I column density, low metallicity DLAs. 

The main problem with both these explanations is that 
they do not explain the striking symmetry that exists 
between the upper and lower envelope, and instead at- 
tribute the symmetry to coincidence. Two possible ex- 
planations that instead take the symmetry as a premise 
are (1) the fundamental plane equation is H I column 
density dependent. If indeed the tilt of the plane is 
slightly different for higher column density systems (i.e. 
their metallicity does not evolve with redshift as much 
and/or their mass- metallicity correlation has a slightly 
different slope) , then this could reduce the scatter in their 
metallicity. Our current data set is lacking enough high 
H I column density DLAs to confirm or refute this expla- 
nation. (2) Higher H I column density systems contain 
more absorption components than lower H I column den- 
sity systems. The thought behind this is that each dark 
matter halo consists of components with varying metal- 
licities. A low H I column density sightline through this 
halo will only sample a few components, and therefore 
can experience a wide range in metallicities. On the other 
hand, a high H I column density sightline will sample 
many components, and its metallicity will be the aver- 
age metallicity of all these components. As a result the 
scatter in metallicities for these large H I column density 
systems is smaller than for the low H I column density 
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Figure 9. Plot of the two kinematic properties where the sample 
has been apportioned into two subgroups based on the redshift 
of the DLA. The low redshift subgroup, solid (purple) squares, 
shows a small scatter of 0.10 dex around a linear fit to the data, 
whereas the higher redshift sample, solid (green) triangles, shows 
a significantly larger scatter of 0.20 dex. This increase in scatter 
is only toward higher Avgo values for a given W1526 value. In the 
text we provide two possible explanations for this observation. 

systems. By varying the average metallicities of each 
dark matter halo depending on redshift and mass, we 
can still reproduce all other correlations with this expla- 
nation. Both explanations are therefore fully consistent 
with current observations. 

Besides the dependency between metallicity and H I 
column density, the second set of notable dependencies 
discussed in this paper are those between redshift and 
the kinematic parameters, Avgo and Wi526- Both kine- 
matic parameters are believed to trace the mass of the 
dark matter halo th at hosts the DLA, although with sig- 
nifica nt scatter (e.g. lLedoux et a l. 2006; Prochask a et al.l 
[20081 see further Section I5.2[) . Therefore, these two 
parameters should show a strong correlation as is dis- 
played in Figure [9l In this figure we have apportioned 
the sample into a low redshift and high redshift sub- 
sample. Interestingly, the high redshift sample has twice 
the RMS scatter compared to the low redshift subsam- 
ple. Moreover, the scatter at high redshift only in- 
creases upward of the AV90-W1526 trendline at low red- 
shift, suggesting that for a given IU1526 value, the max- 
imum available velocity width, Avgo imax , increases with 
increasing redshift. Since the scatter around the W1526- 
metallicity correlation is smaller than the scatter a round 
the Av9o-metallicity correlation. iPro chaska et al. (2008) 
argues that W1526 is a better indicator of the dark matter 
halo mass than Avgo. Under this assumption, the above 
result suggests that for a specific dark matter halo mass 
(i.e. a specific value of Wi526)i Avgo, max increases with 
redshift. This result, however, cannot be statistically 
confirmed with the current sample size. 

Nevertheless, we provide here two potential explana- 
tions for the increase in scatter with increasing redshift 
around the Avgo-Wi526 correlation depending on which 
model is used to describe DLAs. If we assume that 
DLA s are large massive disks (e.g IProc haska & Wolfe 
Il997f ). then one way to interpret this result is that for 
a given dark matter halo mass, the rotational velocity 
structure of the neutral gas increases with increasing red- 
shift, which will increase Avg Uim ax- This is in agreement 
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with the theoretical model of IMo et al.l (|1998l ). which 
predicts that for a given dark matter halo mass the ro- 
tational velocity of the galactic disk will increase with 
increasing redshift. If instead we assume that DLAs 
are a collection of l ow-mass protogalactic clumps (e.g 
lHaehnelt et al.l fl998T ) . as in current ACDM models, we 
can interpret the result as a redistribution of the neu- 
tral gas content in dark matter halos with redshift. In 
this case, the Avgo, m ax for a specific dark matter halo 
mass at high redshift is larger, because some sightlines 
will probe satellite galaxies which move in the poten- 
tial of the primary galaxy increasing the relative velocity 
of these objects. At low redshift, more satellite galax- 
ies will have merged, decreasing the probability of these 
sightlines, and therefore decreasing Avg , ma x- Together 
with a dec rease in galactic scale winds at low redshift, 
Ccn (2012) has shown with numerical simulations that 
these two effects can decrease the Avgo )m ax with decreas- 
ing redshift. 

Although we see tentative evidence for an evolution 
in the scatter around the correlation between Avgo and 
W1526 with redshift, no correlation is found between red- 
shift and Avgo as shown in Figure [5^,- Interestingly, this 
constancy in Avgo wit h redshift i s also seen in the nu- 
merical simulations bv lCenl (|2012f) . iCenl ((2012) attributes 
this constancy to two countering processes: the growth 
of the dark matter halo mass with decreasing redshift, 
and the decrease in galactic scale winds due to reduced 
star formation at lower redshifts. Unlike Avgo, W1526 
shows a potential correlation with redshift. However, 
this correlation cannot be confirmed above the 3-cr level 
with the current sample size. Using a variety of different 
tests, we also looked at any potential evolution in the 
distributions of the kinematical parameters. Both dis- 
tributions are consistent with no evolution. The lack of 
evolution of both the mean and distribution for both pa- 
rameters is significant for two reasons. First, the lack of 
evolution in the distribution of either kinematic param- 
eter indicates that at each redshift, DLAs are embedded 
in a wide range of dark matter halo masses. Second, 
the weak correlation between redshift and kinematics al- 
lows us to combine two separate correlations, namely the 
rcdshift-metallicity correlation and the mass-metallicity 
correlation, into a single fundamental plane equation. 

5.2. Implications of the fundamental plane of DLAs 

The main result in this paper is the existence of a fun- 
damental plane for DLAs. In Section [4] we show that the 
rcdshift-metallicity and metallicity-velocity width cor- 
relations can be combined into a single planar equa- 
tion, [M/H] = (-1.9±0.5) + (0.74±0.21)-logAvg - 
(0.32±0.06)-z. If instead we use the IU1526 statis- 
tic, the equation becomes [M/H] = (-0.47±0.14) + 
(Ll±0.3)-logWi526 - (0.16±0.06)-z. The underlying rea- 
son for the existence of both planes is thought to be 
similar. Since both kinematic parameters are tracers 
of mass, the fundamental plane simply combines the 
rcdshift-metallicity correlation and the mass-metallicity 
correlation, which are correlations known to exist for 
DLAs, into a single equation. 

This equation provides a better description of the fun- 
damental relation that exists between redshift, metal- 
licity and kinematics than the two correlations for two 
reasons. First, it reduces the scatter around the correla- 



tions, and therefore provides a more stringent constraint 
for simulations modeling DLAs. The reduction in scat- 
ter is about 20%, giving a reduced scatter of 0.37 dex 
and 0.38 dex around the correlations. This is still signifi- 
cantly larger than the observational uncertainty for each 
measurement, which is about 0.12 dex. However, we do 
not expect the scatter to reduce to the observational un- 
certainty because a single dark matter halo can produce 
a range in metallicities, velocity widths, and equivalent 
widths depending on where the quasar sightline inter- 
sects the dark matter halo. Hence, a certain amount of 
scatter is inherently part of the quasar absorption line 
experiment. To estimate the size of this scatter, we 
tur n to numerica l simu lations such as those presented 
by iPontzen et al.1 (|2008l ). These simulations suggest a 
wide range of uncertainties depending on which model 
is used to describe DLAs, showing that a scatter of 0.38 
dex can be reproduced from just intersecting dark mat- 
ter halos at differing impact parameters. However, most 
models predict a scatter slightly smaller than this. The 
smaller scatter in the IU1526 plane strengthens the idea 
that this parameter tra ces the kinematics of the halo gas 
(jProchaska et al.ll2008[ ). since the halo gas is assumed to 
be more spherically distributed than the neutral ISM, 
and therefore will be less affected by the impact param- 
eter and inclination angle of the quasar sightline. 

The second reason the fundamental plane description is 
a better description than the two individual correlations 
is that its existence confirms that the dark matter ha- 
los hosting DLA s satisfy a ma s s-met allicity correlation 
at each redshift. iLedoux et alj (|2006l ) showed that this 
held true for their sample apportioned into a low red- 
shift and high redshift subsample. Our result expands 
on their result and indicates that the mass-metallicity 
correlation holds over a larger redshift range. Moreover, 
we show that the zero point of the Avgo-metallicty corre- 
lation is evolving with redshift. Figures [7^, and[SJi show 
that both subgroups of DLAs apportioned by redshift fol- 
low a similar mass-metallicity correlation. This can also 
be seen in Figures Ub and [SJd where the more massive 
dark matter halos (i.e. larger Avgo and larger W1526) 
have at each redshift higher metallicities than the less 
massive dark matter halos. These two figures clearly 
show that at each redshift, the majority of the scatter 
in the metallicity-redshift correlation is caused by the 
existence of a mass-metallicity correlation. This result 
is consistent with the existence of mass-metallicity cor- 
relations found at both low and high r e dshift in star- 
forming galaxies (iTremonti et al.1 120041: ISavaglio et al.1 
[20051 lErb et al.l [2006; Maiolino et al. ri2008(). Although 
these studies used stellar mass and not total dynami- 
cal m ass, there is a strong correl ation expected between 
them (jBrinchmann & Ellisl [2000T ) . 

We can use this dataset to explore the evolution of the 
mass-metallicity correlation for DLAs over the redshift 
range z — 2 to z — 5. At redshifts below 0.75, the stel- 
lar mass-metallicity correlation evolves with redshift in a 
mass-independent way; the slope of the correlation stays 
the same, but the zero p oint of the correlation d ecreases 
with increasing redshift (jMoustakas et al.l 120101 ). Inter- 
estingly, this decrease in zero point and constancy of the 
slope is also seen in Figure [7^ for the Avgo-metallicity 
correlation. Moreover, the decrease in zero point for the 
our DLAs is about 0.3 dex per unit redshift, which is 
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Figure 10. Distribution of the kinematic parameter Avgo for both the obs ervational data and t he numerical model by IPontzen et aTl 
(2008). The solid gray distribution is the distribution of the numerical model by Pontzcn et al. (2008). The hashed (red) sample in panel (a) 
is the complete sample, whereas the hash ed (green) sample in panel (b) is just the subset of DLAs that are part of the low-cool population 
of DLAs as defined in (Wolfe et al. 2008). The left panel shows that the complete sample is inconsistent with the numerical model, since 
a two-sided KS test gives a probability less than 1 X 10~ 6 that the two populations are drawn from the same parent distribution. On 
the other hand, the right panel shows that the low-cool subsample of DLAs is consistent with the velocity distribution found from the 
numerical model. Since low-cool DLAs are believed to arise from less massive dark matter halos, this result corroborates the hypothesis 
that numerical simulations are underproducing the occurrences of DLAs that reside in massive dark matter halos. 



comparable to the decrease in the mas s-metallicity rela- 
tions hip seen at redshifts below 0.75 (Moustakas et al. 
12010ft and between 0.7 and 3.5 (jMaiolino et al.1 120081 ) 
for mass-metallicity correlations determined from stel- 
lar masses. This suggests that DLAs are enriched in 
a similar manner as the star-forming galaxies used to 
determine th e stellar mass-meta llicity correlation. In a 
recent paper. lM0ller et al. (2013) conclude that the zero 
point of the mass-metallicity correlation might not evolve 
past a redshift of 2.6. This is in contrast to our result, 
which shows that the zero point is steadily evolving over 
the redshift range z=2 to z=5. Currently the sample 
of unbiased DLAs is too small to distinguish one sce- 
nario from the other. On the other hand, the zero point 
of the Wi526-metallicity correlation is not evolving with 
redshift. This suggests that the Avgo parameter is a 
better indicator of the total stellar mass inside the dark 
matter halo compared to the W1526 parameter. This is 
not surprising since Avgo describes the kinematics of the 
neutral gas, which should be more strongly correlated to 
the stellar mass than the mass of the dark matter halo. 

5.3. Comparison with Models 

We conclude this paper by comparing the fundamental 
plane equation with numerical models which have specif- 
ically modeled DLAs. To date, DLA simulations have 
tried to reproduce some of the dis tributions of param- 
eters such as meta llicity and Avgo (jPontzen et al.1 [20081 : 
iTescari et al.ll2009T) . and the observed two -parameter cor- 
relations (jPontzen et al1l2008l: lCenll2012t) . but this paper 
introduces a more restrictive criterion which should be 
satisfied, the fundamental plane of DLAs. A successful 
simulation should be able to populate the fundamental 
plane in a similar manner to what is observed. How- 
ever, currently no simulation tracks these three parame- 
ters accurately. Instead we therefore consider only pro- 
jections of the fundamental plane. In particular, the pro- 
jection along the Avgo axis gives the metallicity-redshift 
correlation; its agreement with models has been dis- 



cussed in detail in iRafelski et al.l ([2012). The projec- 
tion along the redshift axis leads to the metallicity-Avgo 
correlation. Thi s correlation has bee n co mpared to sim- 
ulatio ns in both ITescari et aTl ()2009l ) and IPontzen et aTl 
(2008), where the latter finds a correlation equation of: 
logi Av 90 = 2.5 + 0.58-[M/H] at redshift 3. They com- 
pared this to the o bservational correlation equation of 
iLedoux et al.l (|2006l) and noted the good agreement ex- 
cept for the significant ly larger obse r vation al scatter. 
Comparing the model of lPontzen et al.l (j2008l ) to the new 
fundamental plane equation at redshift 3 (i.e. setting 
z = 3 in the fundamental plane equation) gives a sim- 
ilar agreement. The main improvement here, however, 
is in the scatter around this relationship. The scatter in 
previous observations was too large, because these obser- 
vational analyses did not take into account the redshift 
evolution of metallicity. Using the fundamental plane re- 
duces the scatter by abou t 20% , and therefore the scat- 
ter in the IPontzen et al.l (j2008l ) simulation is in better 
agreement with the observational data then previously 
assumed. 

Finally, we also compare the distributions of each 
of the parameters. Of these, the Av 90 distribution 
is the most interesting, because of the long stand- 
ing i nability of models i n rep r oducing this di s tribu- 
tion (iProchaska fc Wolfd Il997t iHaehnelt et~aTI 119981: 
iRazoumov et al.l 120061 : IPontzen et al.ll2008| ). As was dis- 
cussed in the introduction and shown in Section 13.11 
this distribution has two characteristics that are difficult 
to reproduce by simulations based on ACDM models; 
namely the large median velocity, and the high- velocity 
tail of the distribution. It is believed that the inability 
of reproducing these features of the Avgo distribution 
is due to the inability of the numerical simulations to 
produce enough DLAs that reside in massive dark mat- 
ter halos compared to less massive dar k matter halos 
(jPontzen et al l 120081: ITescari et all I2009D . We test this 
hypothesis by comparing the numerical simulations to 
just those DLAs that are believed to reside in smaller 
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dark matter halos. 

To do this we select only the PL As from our sa mple 
that have low cooling rates, l c . I Wolfe etldl ()2008l ) sug- 
gests that these DLAs reside in less massive dark matter 
halos compared to the high cooling rate DLAs. Figure 
[TUb shows the Avgo distribution of the complete sam- 
ple (shaded red histogram); the gray solid histogram 
is the distribution obtained from the simulations by 
iPontzen et al.l (|2008[ ) . As can be seen from the inset and 
the KS test probability, there is a less than 1 x 10~ 6 prob- 
ability the two distributions are drawn from the same 
parent population. On the other hand, Figure [TUb shows 
that the DLAs with a low cooling rate (shaded green 
histogram) have a Avgo distribution similar to that ob- 
tained from the numerical simulation. This result sup- 
ports the hypothesis that the numerical simulations are 
capable of reproducing the Avgo distribution of the DLAs 
in less massive dark matter halos, but have difficulty pro- 
ducing enough DLAs that reside in more massive dark 
m atter halos. This conclus ion was also recently drawn 
in iFont-Ribera et al] (|2012f) using the cross-correlation 
between DLAs and the Lyman a forest. 

One way numerical simulations try to mitigate the dis- 
crepancy in the observed and simulated velocity width 
distribution is by adding in galactic scale outflows 
(jTescari et all 12003 : lHongeTal1l2(n"ol lc!e^ l2CTT2h . These 
outflows are capable of decreasing the cross section of 
less massive dark matter halos, somewhat decoupling the 
Avgo statistic from the mass, and increasing the cross 
section of the more massive dark matter halos. All of 
these effects increase the relative number of DLAs in 
more massive dark matter halos. With the addition of 
galactic scale outflows, these models have been relatively 
successful in reproducing the Aygp distribution o f DLAs 
(|Tescari et al.ll200l lHonglFaIll2010l: ICe^l[20T2h . which 
is supported by the sa mple in this paper. However, 
IFont-Ribera et al] (|2012l ) suggest that this process still 
might not produce enough DLAs hosted within massive 
dark matter halos, indicating that there is still some dis- 
crepancy between the observations and the simulations. 

6. SUMMARY 

We have studied the spectra of 100 DLAs taken with 
a single high resolution instrument, the High Resolution 
Echelle Spectrometer on the Keck telescope, to deter- 
mine the fundamental relations that exist between the 
measured parameters (Table [4]). We find the following 
new results. 

1. The metallicity-H I dependency shown in Fig- 
ure [3] shows a lack of high H I column density, 
low metallicity systems besides the known lack of 
high H I column d ensity, high metallicity systems 
(Boisse et al. 1998). Possible explanations for the 
lack of high H I column density, low metallicity 
systems are that (1) low metallcity systems arise 
from cold flows which are not dens e enough to form 
high H I column density systems (jFumagalli et al] 
120111 ). or (2) that higher H I column density sys- 
tems arise from sightlines probing the inner part of 
galaxies, which have higher me tallicity because o f 
an existing metallicity gradient (jChen et al.ll2005l ). 
However, neither of these explanations can explain 
the reflection symmetry that exists about the line 



[M/H] = —1.43. One possible explanation that 
takes this reflection symmetry as a premise, is 
that higher H I column density systems probe a 
larger number of components. The larger number 
of components results in a smaller probability of 
encountering a sightline with a very low or very 
high metallicity. Hence, this will decrease the vari- 
ance around the mean metallicity value of —1.43 
for higher column density systems. 

2. The Avgo parameter is not evolving between red- 
shift 2 through 5 (Figure [5b) • This indicates that 
at all redshifts, DLAs are hosted by a wide vari- 
ety of dark matter halo masses, and it allows for a 
combined description of the Avgo-metallicity and 
redshift-metallicity correlations. 

3. The main result of the paper is that we can de- 
scribe the relations that exist between the redshift, 
metallicity and mass of a DLA by a single fun- 
damental plane equation: [M/H] = (— 1.9±0.5) + 
(0.74±0.21)-logAv 90 - (0.32±0.06)-z. This plane 
equation has as advantage that it reduces the scat- 
ter around either of the correlations; providing 
a more stringent constraint for numerical simu- 
lations. Secondly, it confirms the existence of a 
mass-metallicity relationship at each redshift be- 
tween redshift 2 through 5, where the zero point 
evolves with redshift. This evolution in the zero 
point with redshift is consistent with the evolution 
of the zero point of the mass-metallicity relation- 
ship seen in star forming galaxies. 

4. Finally, we compare the sample data to numerical 
models, and find that numerical models are unable 
to reproduce the Av 9 o distribution of the complete 
sample (Figure [TUb). However, these models are 
able to reproduce the Avgo distribution of the low- 
cool subset of the sample (Figure \T0b) . These low- 
cool DLAs a re believed to res ide in smaller dark 
matter halos (| Wolfe et al.ll2008f > . This result there- 
fore supports the hypothesis that numerical models 
fail to pro duce enough DLAs i n massive dark mat - 
ter halos (|Pontzen et alj|2008t pTescari et al.ll2009D . 
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